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1. Introduction 

A new and exciting field in condensed matter physics started when graphene - a two- 
dimensional, one carbon-atom thick material - was isolated for the first time.|[Il 13 It was 
experimentally shown that the charge carriers in graphene could be controlled by a bottom gate 
setup-up; the charge carrier were shown to be either holes or electrons depending on the sign 
of the bottom gate voltage. In the transition from hole-based to electron-based transport the 
conductivity shows a minimum (not zero) value, o"min. Its experimental value is of the order 
of o"min — 4e^//i,|[Ill2l[3l|4l but the actual value seems to be some what sample dependent. jH 
This value for cTmin imposes therefore a limitation on the minimum value of the current a field 
effect transistor made of graphene can transport. The existence of a conductivity minimum in 
graphene is a consequence of the fact that the elementary excitations of graphene are Dirac 
fermions, with a linear dispersion relation, instead of usual electrons with parabolic-like 
dispersion, characteristic of ordinary semi-conductors. Interestingly enough, the calculated 
value of the conductivity of graphene at the neutrality point is off the experimental value by 
the factor l/vr. [8l| Although this value is the more common result for the theoretical 

calculation of a^i^, there are, however, several different values available in the literature. Il9l 
It is also interesting that a clean graphene sample with metallic leads and smooth edges has 
a value of cr = al^^° as long as it width (w) is much larger than its length (L), being smaller 
than a^^° in the opposite limit. ifTOl Considering the case of metallic armchair edges, it found 
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that a > a^^^° for w/L < 1 and that a cr^™ for w/L:^ l.iM This shows that disorder 
is not needed for having a ~ o-^^in- 

Another characteristic of Dirac electrons in graphene is their ability to tunnel through 
a potential barrier with probability one. [[TTl \12] This so called Klein tunneling of chiral 
particles has long ago been proposed in the framework of quantum electrodynamics . [1131 [T4l 
[TSll but was never observed experimentally. Graphene opens up a route to observe this effect 
in a tabletop experiment, where the potential is created by some electrostatic gate potential. 
The manifestation of Klein tunneling is also present when electrons in graphene are forced to 
transverse a n — p junction, leading to a selective transmission of those electrons approaching 
the junction perpendicularly. [16] Other unusual effects, such as the focusing of an electric 
current by a single p — n junction are also characteristic of Dirac electrons in graphene. ifTTl 

As appealing as the Klein tunneling may sound from the point of view of fundamental 
research, its presence in graphene is unwanted when it comes to applications of graphene 
to nanoelectronics. This comes about because the pinch-off the field effect transistor may be 
very ineffective. The same may occur because of the minimum conductivity of graphene at the 
neutrality point (as discussed above). One way to overcome these difficulties is by generating 
a gap in the spectrum. From the point of view of Dirac fermions this is equivalent to the 
generation of a mass term. There are two known forms of generating gaps in the spectrum 
of graphene. The first one is by patterning graphene nanoribbons. niSl iT9l The mechanism 
of producing these gaps depends on the nature of the termination of these nanoribbons. For 
armchair nanoribbons the gap comes from quantum confinement of Dirac fermions induced 
by the finite nature of the ribbons in the transverse direction. For zigzag nanoribbons the gap 
stems from the formation of polarized spin edge- states characteristic of these type of ribbons. 
The formation of these polarized states is also possible in bilayer graphene. [20] It is interesting 
to notice that Klein tunneling can also be circunvented by using a graphene bilayer. |[TT]| The 
value of the induced gaps depend on the width of the ribbons, but for large widths it is of the 
order of 0. 1 e V. 

Another possibility of generating gaps in the graphene spectrum is to deposit graphene 
on top of hexagonal boron nitride (BN). [|21 | | This material is a band gap insulator with a boron 
to nitrogen distance of the order of 1 .45 A, [[23l (in graphene the carbon-carbon distance is 
1 .42 A) and a gap of the order of 4 e V. It was shown that in the most stable configuration, 
where a carbon is on top of a boron and the other carbon in the unit cell is centered above a 
BN ring, the value of the induced gap is of the order of 53 meV. Depositing graphene on a 
metal surface with a BN buffer layer leads to n— doped graphene with an energy gap of 0.5 
eV.[|22l 

The two mechanisms described above can be used to produce arrangements of graphene 
where in some spatial zones of the material the Dirac electrons will have gaps in the 
spectrum. The first possibility is to pattern graphene planes such that in several areas of 
the graphene flake narrow nanoribbons may exist. Another possibility it to combine wafers 
of silicon oxide and hexagonal boron nitride, such that in the region where the BN is 
located the local spectrum of graphene will present a finite gap. We shall explore in this 
paper this latter possibility and the way it can prevent Klein tunneling from occurring. The 
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two mechanisms just mentioned can then be at the heart of future nanoelectronics built of 
graphene. The second method is also related to junctions of graphene with other kind of 
systems, being them superconducting. [ [24l l25l normal-conductor/graphene/normal-conductor, 
ll26l or multi terminal j unctions . [jlTI Also the study of electron transport in disordered 
graphene samples is of interest, [281 specially because the tunneling may be assisted by 
impurities, [|29l which is a manifestation of Klein tunneling. 

For a review on the experimental aspects of graphene physics see the work of Geim and 
Novoselov. ll30ll Some of the theoretical aspects of graphene physics are reviewed qualitatively 
by Castro Neto et aZ.. [[311 by Katsnelson, [[321 and by Geim and MacDonald; [33] a more 
comprehensive review is given by Castro Neto et a/.. [[341 For a review on Klein tunneling see 
the work by Beenakker. [[35l 

2. Basic definitions. 

As described in the previous section, we assume that it is possible to manufacture slabs with 
Si02-BN interfaces, on top of which a graphene flake is deposit. This will induce spatial 
regions where graphene has a vanishing gap intercalated with regions where the BN will 
cause a finite gap. 

In the following we will consider the graphene physics in two different regions: the 
A;— region, where the graphene sheet is standing on top of Si02, and a g— region, where a 
mass-like term is present, caused by BN, inducing an energy gap of value 2t' (for all numerical 
purposes we use t' =0.1 eV). The wavefunctions in these two regions will be referred by 
and ipq, respectively. The geometry of the scattering process is represented in Fig. [T] 

The Hamiltonian for massless Dirac electrons in graphene, around the 1^ -point in the 
Brillouin zone is given by 



where cr = (cr^, (7y),p = —ih'V, ai, with i = x,y, z, is the i Pauli matrix, and vp = 3ta/{2h), 
with t the nearest neighbor hopping matrix in graphene and a the carbon-carbon distance. 
Therefore, in the massless wave function, in the A; —region (t' = 0), is given by 



Hg = vpo- ■ p 



(1) 




(2) 



with 



M(k, s) = se , 

s = sign(£') and 9 = arctan (ky/k^). The corresponding energy eingenvalue is 



(3) 




(4) 



with k the absolute value of the wavevector. 

In a region of finite mass the Hamiltonian for Dirac electrons is 

Hg = vpCT ■ p + t'a^ , 



(5) 
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with mv'^p = t' the mass term (m is the effective mass); as a consequence the electronic 
spectrum will present a finite energy gap of value 2t' . In the g-region (the gaped region, 
t' 7^ 0), the wave funcion is 

= ^ f / e^''\ (6) 

where 

^/ 

^(q, s) = - — — . (7) 

nvpigx - iqy) 

Due to momentum conservation, electrons propagating through ak — q interface will conserve 
their wavevector component parallel to the interface. Thus, taking this interface to be located 
along the y axis, we will have always ky = qy. The g-region eigenenergy, associated with the 
eigenstate ^ and the Hamiltonian ([5]), is 



E = ±^{q^^ + kl){hvpy + t'^. (8) 

It is amusing to notice that the spectrum ([8]) has the same form as for the electrons in 

a graphene bilayer, when the two graphene planes are at different electrostatic potentials. 
[[361 [371 Using Eqs. ^ and ©, we write 



vphq^ = cos2 (6) - t'2 (9) 

and, depending on E'^ cos^ (9) being larger or smaller than t'"^, q^ may take a real or a pure 
imaginary value. Wave propagation follows for the former case, evanescent waves in the 
latter. 

For a real g^, and since qy = ky, we have 



_ t'2 sin (0) = 1^1 sin {6) (10) 

where is the angle of propagation of the electron in the g-region (see Fig[T])- Equation ([TO] ) 
is just the usual Snell's law, for electrons being refracted at the interface separating the k— 
and g— regions. We see that 4>>6 whenever \E\ > t' . 



2.1. Forward and backward propagation. 

We consider now the simple reflection in the interface, with the incident and the reflected 
waves both on the k— or on the g— region. 

In the /c— region, the x— component of the wavevector of the reflected wave is 
symmetrical with respect to the incident wave. Thus, for this case, we have the following 
transformations under a reflection (see also Fig{T]) 

k^^-k^ and e^^ e'("-^) = -e-'^ (11) 

This leads to the generalization of Eqs. Q and (|3]) 

M± = ±se±*^ (12) 



(13) 
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Figure 1. Top figure: graphene band structure for the massless and massive cases. In this latter, 
the quasi-parabohc bands have a gap energy 2t' . Bottom figure: geometry in the reflection 
in the k — q interface. An incident wavefunction ?/;^ with wavevector k+ is reflected and 
refracted into the ip'j^ and the ?/;+ wavefunctions with wavevectors k_ and q+, respectively. 
Since the momentum is conserved at the interface, one has that q^, = kj,. The refracted wave 
propagates with an angle cf), which is slightly larger than the incident and reflected angles 9 
with Iq^:! < |ka;|, a consequence of energy conservation. 



where the plus and minus signs refers to waves propagating, respectively, in the positive and 
negative directions of the x-axis. 

A similar reasoning leads to the generalization of ^ to v{q, s), 

^± = t. H"^' -i-v ^^^^ 
nvF[±qx - iky) 

and, also, of the g— region wavefunction to 

The differences we have just highlighted on the wavefunctions and coefficients for 
forward and backward propagating particles can be also seen in the differences in positive 
and negative angles of incidence in the interface. This changes are useful, when a guiding- 
wave kind of device is made. Let us therefore analyze the case when ky —ky. If in Eq. 
(fTTI) we keep kj. unaltered and "reflect" instead ky we would obtain 

ky -ky and e^^ e"'^ (16) 

with similar relations for (p, the angle in the q-region. For this cases, we get 

v±{-<p) = ) and n±{-e) = -n^{e). (17) 

Apart from a minus sign, these relations shows that the the operation of changing the sign of 
ky (i.e., the angle of incidence) is equivalent to the one of changing the sign of k^ (q^ in the 
q-region). Of curse that the extra minus sign in the right hand side of both expressions in Eq. 
(flTl) are of no consequence within the calculation of reflection and transmission coefficients 
that follow. 
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2.2. Real and evanescent waves in the q-region 

Since in the g-region there is a gap in the energy spectrum then can take both real and pure 
imaginary values. In the first case, we have wave propagation in this region, in the latter just 
evanescent waves. No simple expression as the one given by Eq. (fT2)) can be written in this 
case. Instead, we need to consider separately the cases where q^ is a real or a pure imaginary 
number. 



2.2.1. For q^ real. For real q^, we can write a similar expression to the one in Eq. (fT2l) . 

v± = ±ve^'^, (18) 
with 4> given by Eq. (fTOl) and 



_E-t' _ E-t' 
where Eq. ([8]) was used. 

2.2.2. For q^ pure imaginary. Since ky is always a real number, Eq. (fT4l) implies that, if q^ 
is a pure-imaginary, v± also is, and then 

v± = TiJ^±, (20) 

where 

^± = ± ^ f,"/ V (21) 
nvpip^Ky — Oi) 

with the real absorption coefficient a defined as = za, and a given by 



a = {vphy^^t'^ - ^2 cos(^)2 . (22) 

2.2.5. Complex conjugate of the u± and v± coefficients. For the calculation of the intensity 
reflection and transmission coefficient we will need to deal with the complex conjugate of the 
u± and v± coefficients. The definition (fT2l) for u± implies that 

u*j^ = — n=p. (23) 

In the case of v±, its complex conjugate depends on the the fact of having a real or 
imaginary q^, 

f —v^ if qx is real 
v± = { ^ . . . (24) 

[ —v± it qx IS imaginary 

3. Transmission and reflection at the interface: the step case. 

3.1. Refiection and transmission amplitude coefficients. 

We compute now the refiection and transmission amplitude coefficients for electrons crossing 
an interface between a k— and a g— region. Unlike what happens in optics and due to the 
differences on back and forward propagation, we will need to consider not only two but four 
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Figure 2. The four different possible cases for reflection/transmission in an interface between 
k and q regions. 

different cases: electrons crossing the interface coming from the fc— region in the forward and 
backward senses, and those crossing the interface coming from the g— region, also propagating 
on the positive and negative senses of the x axis. These four cases are summarized in Fig. [21 



3.1.1. Propagation from a k— into a q— region. We start by deriving the amplitude reflection 
and transmission coefficients, which will be denoted as r^^ and t^^ respectively, for the case 
of the propagation from a k— into a q— region. This situation is described in Fig. [T]and also 
in Fig. |2la). 

Since there is a partially reflected wave, the total wave function in the A;— region must be 
written as an superposition of one associated with the incident electrons and other with those 
that are reflected, 

vl/fc(r) =A^+ + i?V'k- (25) 

A and B are the normalized amplitudes for the incident and reflected wave functions. In the 
g-region, with C the amplitude of the transmitted wave function, we have 

^,(r) = C^+. (26) 

Using Eqs. (l25l) and (|26l) . and imposing the continuity condition of the particle's wave 
function at the interface, i.e. ^ki^, y = 0) = '^q{x, y = 0), we find 
B _ v+~u+ C 
A- u^-v^ ^"""^ ^'^^"A 
where the superscript + recall that the incident wave function is, in this case, traveling in the 
positive direction of the x-axis. 

Had we considered the case were the particles travel in the backward direction, 
represented by Figl2lc), we would have obtained 

_ V- — U- 

^kq = ^"^^ ^fcg = 1 + ^kq- (28) 

This result can be obtain simply by exchanging the plus by the minus signs in Eq. (1271 ). 



rtq = ^ = — and tj^^ = - = 1 + ij^^, (27) 
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3.1.2. Propagation from the g— into the k— region. For computing the reflection and 
transmission coefficients for the cases where the electrons come from the g— region into the 
fc-region, r^, and t^^,, we need only to exchange u ^ v m the corresponding backward and 
forward expressions (l27l) and (|28l) . The result is 



3.2. Amplitude coefficients: general algebraic relations. 

It will be very useful in the following, for expression simplification purposes, the derivation 
of simple relations between the reflection and transmission amplitude coefficients. Similar 
relations to those we present here exist also for the photons' optics case. For instance, we 
may write ±ri2 + ti2 = 1 when a light beam is reflected and refracted in a diopter between 
regions 1 and 2, with the plus or minus sign corresponding respectively to the cases where ni, 
the index of refraction of medium 1, is smaller or larger that the one of region 2. 

Here, however, we have that in general r+„ 7^ (and similarly for the transmission 
coefficients) and these relations are less trivial (we have used the notation m 7^ n = {A;, q]). 

Using the definitions in Eqs. (ITTI) . (l28l) and (|29l) . we can write 

7^ + T = l, (30) 

where 

' <Y>~^ <Y> <Y>~^ ly 

— ' kq' kq — ' qk' qk /oix 

r -f+f+ - f-f- ^ ^ 

^ ''qk^kq ^qk^kq- 

These relations are general and don't depend on the being real or imaginary. Another 
general relations, useful to simplify expressions of the transmission of multi-layered 
structures, are 

^± _ _7? XX rim, /o^x 

' mn' nm f± ' ^ 

nm 

3.3. Intensity reflection and transmission coefficients. 

The general definitions for the intensity reflection and transmission coefficients are 

(33) 



mn \' mnJ 



mn mn \ mnJ ''mni 

where we keep the same notation as before. We will consider now, separately, the cases were 
is a real number or a pure imaginary. 

3.3.1. For q^ real. For q^ a real number, we note first that for any m 7^ n = {/c, g}, 

(34) 



V mnl mn 



Tz = n* 



Tunneling ofDirac electrons through spatial regions of finite mass 9 

These relations are just a consequence of the definitions of in Eqs. (|27l) . (|28] ). and (|29l) . 
and in Eqs. dlS]) and (EH). Using Eq. dH) in Eq. ([33]) results in 

R = TZ and T = T for g^, real, 

which is valid, for both interfaces and both directions of propagation. Furthermore, using Eq. 
(I30l) we get 

R + T = l, 

an expected result. 

Explicitly, the TZ coefficient defined in Eq. (|3T1) is given by 



^ ^ {v+v- - 1) - 



Making use of Eqs. (|T4l) and (|T2l) . we may write 

and u±v± = sve^''-^+'^\ 
where v is given by Eq. (|7]) and (p by Eq. (flOl ). Using these expressions we obtain 

1 + f ^ — 2 s f cos (6' — 6) 



n 



1 + + 2 s u cos (6' + 0) ' 
where 

w cos (6* ± 0) = (gx' cos (6') =f ky sin (6*)) . 

h/ -\- 1 

Finally, after algebraic simplification, we obtain 

R = n= (35) 

kx + Qx 

3.3.2. For Qx pure imaginary. For a pure imaginary, we see that 

fr^ = 1 

\ mn/ mn 

(36) 

nn* = 1 

Using these relations along with Eq. (|33l) . we straightforwardly obtain 

= 1 and T^g = with Qx a pure imaginary. (37) 

This is an expected result since the transmission T = 1 — R must be zero in the case where 
the wave in the A; -region enters in the gap of the g-region. If the incident wave propagates 
in the gap region, i.e. it is an evanescent wave, the coefficients Rff. and are physically 
meaningless. 

We see from the second expression in (|36l) that 7^ is a modulo 1 complex quantity. It 
may be written as 

n = e''^'^ , (38) 
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with 2^9 a convenient definition of its argument. To compute this angle, in the spirit of Eq. 
(|22|) . we replace by ia in Eq. (l35l) to obtain 

kx — ia 



n 



kx + ia 



Computing the real part of this quantity, we get 



t'2 



COS {2if) 

and, after straightforward manipulation, 

cos {if) t' = cos (9) \E\ or else tan 



a 
kx 



(39) 



a Snell type expression for the qx pure imaginary case. 

Since R = TZTZ* = 1 for qx a pure imaginary, a general expression for the intensity 
reflection and transmissions coefficients (valid for g^. both real and pure imaginary) is given 
by 

i?= 1-T 

kx qx 



kx + qx 



2s f cos (9 — 



l + i;2 + 2sf cos(0+ ' ^^^'^ 
Naturally, Eq. (l40l) depends on t', since both v, 9, and depends on this quantity. When 
one considers the case 9 = (f) = t' = one obtains R = 0. This expression is plotted in a 
density /contour plot in FiglH 



4. The barrier. 



With the above definitions, the computation of transmission and reflections coefficients for 
any type of multi-interface device follows similar expressions as those found in normal 
opticslj To illustrate this, we consider in the following a hetero structure made of a g— region 
of width w placed between two semi-infinite slabs of fc-regions, as shown in the Fig. HI Our 
goal will be the derivation of the intensity transmission coefficient for this case, which we will 
denoted by T^. We notice that results for barriers of the same height when the spectrum of the 
the electrons is linear in every spatial regions was considered in Ref. f38l. 

In FigureHl the wave function describes an electron, traveling in the positive direction 
of X— axis, just before crossing the diopter q — k. This wave function can be seen as resulting 
from the coherent superposition of two wave functions, one being itself after a round trip in 
the g— region, given by tjfct"^ ^^'^qxw^ another one which is the incident wave funtion 
^0 after crossing the first interface k — q, equal to t^q e*'^^"'. Adding these two contribution 
and solving in order to we obtain, 

= \l>o '"^ — ■■ • 

' qk' gk^ 

I There is no analog, however, for the gap-region with normal incidence. 
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Figure 3. Intensity transmission for particles crossing the interface from a fc-region into a 
g-region. The black region corresponds to a zero transmission, a case that correspond to the 
total internal reflection in the usual photonic optics. 
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Figure 4. Barrier: scheme for the computation of the transmission. 



If we denote the amplitude transmission coefficient for this barrier as 
the fact that \l/2 = ^^fc^i^ we finally obtain 

tn = 

1 - 7^e*29^"' 

where the definitions (l30l) were used. 



^2/^0 and using 
(41) 



4.1. qx real: free propagation. 

If there is wave propagation in the g-region, is real, IZ = R and T = T, and 



tRt%, = 1 — R 



B — i^B i^B 



B 



1 + ( ^J^^ sin^ {q^w) 



where we used the finesse definition 

^/R _nt' 



J^ = vr- 



(42) 



(43) 
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Figure 5. Transmission of a barrier for a q-region width of w — 50a (top) and w = 300a 
(bottom). For a sufficient narrow width, the wave tunnels across the g-region resuhing in a 
non null transmission. In optics, this behavior is known ?ls frustrated total internal reflection. 
The dashed lines marks the region where in the step, the transmission is zero. 

to highlight the similarity with an Fabry-Perot solid etalon (made of glass, for example) in 
usual optics. [39 1 However, this similarity is elusive. In the solid etalon case, in general, the 
finesse is almost a constant coefficient since the interfaces' reflectivities (e.g., in a diopter 
glass-air) has a small dependence on the energy (optical resonances are typically far from the 
visible part of the spectrum and there's no gap as in the case of the graphene with a mass 
term). In the case treated here, JF has a strong dependence on the energy E of the particles 
and, furthermore, there is also a gap present. We will revisit a Fabry-Perot type of device 
further in this work. 

4.2. Inside the gap: frustrated total internal reflection. 

Inside the gap, is a pure imaginary and there's no wave propagation. This is similar to the 
total internal reflection in optics, where only an evanescent wave exists that carries no energy 
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Figure 6. Transmission coefficient for a simple step and several barriers for normal incidence. 
In the barriers case, there is a frustrated total internal reflection and, in the gap, the 
transmission is non zero and increases with decreasing values of the width of the q-region. 




-10 -5 5 10 

E/f 



Figure 7. Left: Transmission coefficient for a barrier with a width 50a and for different angles 
of incidence 9. Ai E ~ 0, all curves have the same value and for 6 = i: jl, the transmission is 
a delta function at = 0. Right; Zero energy transmission coefficient of a barrier as function 
of its width w. 

(since in here, the coefficient R^q = 1) and decays exponentially in the x direction (although 
keeping the phase term e*'^" ^ in the y direction). 

However, by placing a /c-region nearby the evanescent wave, some of the energy of the 
totally reflected wave tunnels throughout the gap region, a phenomena known in optics as 
frustrated total internal reflection. This phenomena is also described by Eq. (|42|) . whenever 
qx becomes a pure imaginary. In this case, replacing = ia and using the definitions (|22]) 
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Figure 8. Schemes for the computation of the transfer matrices. 



C = — rTT-v^> (45) 



and (|38l) we may simplify Eq. (|42l) to 

TB = tBt*B = l-RB = -^^, (44) 

where we have used the definition 

sinh^ {aw) 
siv? {if) 

which will be used in the following. If _E = 0, Eq. (l39l) implies that Lp = — 7r/2 and 
a = t'/vph. Tb is in this case independent of 9 and is equal to 

Tb{E = 0) = cosh"^ [w^]. (46) 

This behavior is clearly shown in left panel of Fig. Ul In the right panel of Figj7]it is shown 
how this tunneling transmittance at zero energy varies with the barrier width. A 50% reduction 
is accomplished for a barrier with a width of approximately 36 a. 



5. Transfer matrices 



The method used in the last Section for computing Ib, although being simple becomes very 
difficult to handle for more complex hetero- structures, with more than two interfaces. These 
type of cases are usually treated with the use of transfer matrices. These will be computed in 
the following. Figure [8] shows the scheme used to compute the transfer matrix in an interface 
k — q and q — k.ln both cases, our goal is to derive and from the knowledge of \l/o and 
%. Defining 
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Figure 9. The diode heterostructure: two thin slabs of q— regions of width w separated by a 
fc— region of width d, all inside semi-infinite slabs of fc— regions. 



where Mm„ is the transfer matrix for the generic m — n interface, and using Eq. (I3TI) . we 
obtain the result 



1 



t: 



t: 



(47) 



The determinant of this matrix is given by 



Det{M, 



kq) 



t 



(48) 



As expected, DetCWLkq x Mg^) = 1. The free propagation of a particle in a k— and in a 
g— region of width ^ is given, respectively, by 



















(49) 



6. The diode. 



We consider now a more complex system composed by a sandwich of two g-regions of width 
w separated by a slab of a /c— region with width c?, inside two semi-infinite /c— regions. To 
derive the amplitude transmission coefficient of such a structure we need to compute the 
expression. 

The result of this expression can be simplified using Eq. (|32l) . resulting in 

^ ~ iJZ e^^'i"'^ — ly — 7^(e^*^^"' — l)2e2jA:a;rf' 
For the most important case where the g— regions are barriers with energy higher than the 
energy of the particles, we have a resonant diode. In this case and using the definitions (|22)) 
and (|38l) we get 



sinh 



\aw] 



A2kxd 



sinh^ {aw + ii-p) 



(50) 
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Figure 10. Transmission of a diode structure with w = 50a and d = 100a (top) and d = 200a 
(bottom). 



where is the amplitude transmission for a simple barrier, given by Eq. (|4TI) . 

We may simplify Eq. (l50l) by expanding the term sinh (aw + i^p) and expressing the 
result in a complex polar representation. Doing this and using the definition for C, in Eq. (1451) 
we obtain 



sinh (aw) 



c 



sinh^ {aw + iip) 1 + C 
with the phase term argument given by 



exp 



(p = — arctan [coth {aw) tan (f/?)]. 

The intensity transmission coefficient can now be easily computed, being equal to 

1 



where now, the diode finesse Tn is given by 



^,, = 7rVC(l + C). 



(51) 



(52) 
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6.1. Revisiting the Fabry-Perot: etalon made with "mirrors". 

The expression (|5T1) results in the simple case of a Fabry-Perot etalon if: 

(i) aw ^ 1, which implies coth {aw) = 1 and (p = Lp; 

(ii) E <^ t', which implies that (p ^ 7r/2. 

With these approximations we get 

71 

J-'d = — sinh {2aw) 

and then 

^ 1 + sinh^ {2aw) cos^ (A^^jC?) ' ^ ^ 

an expression that can be derived from a delta function potentials treatment, treated in Section 

m 

6.2. The diode tunneling conductance. 

Let us now compute the tunneling conductance of the device as function of the potential bias 
V , the chemical potential of the leads, and the length w of the barrier. We shall assume that the 
device is operating in the region where the chemical potential of the leads lies inside the gap 
of the barrier. The total tunneling current density (i. e. the current per unit of cross-sectional 
length) through the device is given by 

2e 

X [f{E - fiL) - f{E - fin)] (54) 

where/(x) = + /i^ = fi + eV/2, fiji = yU — e\^/2, and the chemical potential 

of leads in the equilibrium state. The linear response conductance per unit of cross-sectional 
length is given, at zero temperature, by 



J{V,w) = / dOEdE cose T{E, 



^2 



1^1 



tt/2 



Gifi,w) = ^TT^ / decoseTifi,e). (55) 

In Fig. (fTTI) we plot G{^, w) as function of fi for several widths w. It is clear that the value of 
G{fi, w) may change by several orders of magnitude, close to /i = by a small change of /i. 
Naturally for wider barriers one obtains a smaller conductance. 



6.3. A limiting case. 

The limiting case of barrier can be represented by a delta-function potential, V{x,y) = 
gaz5{x). (See next section for complete discussion of delta function potentials in the Dirac 
equation.) It is interesting to compute the reflected flux for both Schodinger and Dirac 
electrons for this potential. In the first case one obtains 
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10 



-2.0 



-1.0 



0.0 

^l/t' 



Figure 11. (color on line) Linear response conductance G(/i, w), per Dirac cone, as function 
of the chemical potential ji, for different values of the width w. The hopping matrix t is taken 
to be t = 2.7 eV and t' 0. 1 eV. 



whereas for Dirac electrons the result is 



(57) 



It is clear that for electrons in graphene R is angular and energy independent. For the case 
g ^ hvp the reflexion tends to unity. 

7. The diode: a limiting case. 

Finally we want to discuss a limiting case of the resonant tunneling diode made of graphene. 
The device is represented in Fig. [9l The corresponding study for Schrodinger electrons was 
done by Tsu and Esaki. [|43l 

A limiting situation of the device described in Fig. |9]is one where the barriers are are 
describe by a scalar Lorentz potential of the form 



The connection with the true barrier is made by identifying g with at'wa, with a a numerical 
constant of dimensions inverse of length. This form of the potential is equivalent to a mass 
term and therefore to a gap in the spectrum. However, given the short range nature of the 
potential, its effect comes only in the boundary conditions imposed on the wave function at 
the potential position. 

The problem of Dirac electrons in delta function potentials has been studied in the 
pastdm HH and is not without subtleties. Il46l l47l l48l The subtleties can be traced back 
to the problem of evaluating the integral 



V{x, y) = hm^ — [1 - 9{\x\ - e)](T; 



+ ]irng-[l-9{\x-d\-e)]a 



= gaz[S{x) + 6{x — d)] . 



(58) 




(59) 
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where /(x) is a discontinuous function at x = 0. If we try to solve the problem of Dirac 
electrons with a delta function potential using the same trick ll42l one uses for Schrodinger 
electrons we face the problem defined by the integral (l59l) . This is so because the wave 
function of Dirac electrons in a delta function potential is discontinuous at the point where 
the delta function is located. There are several strategies to overcome this difficulty. 
Il44l l45l l46l l47l l48l The most straightforward was devised by McKellar and Stephenson 
||44l|45]| and generalized by Dominguez-Adame and Macia.03 In short, the Dirac equation 
along the x direction can be can be written as 

^ = , (60) 

where (f){x) is a spinor wave function. This problem can be formally solved as 

0(x) = T.e^^'o ^^^^'^Xxo) , (61) 
where the operator is the position order operator such that 

T,[Gix)Giy)] = Gix)Giy)eix - y) + - x) . (62) 

Since we are interested in determine the boundary conditions obeyed by the wave function 
(f){x) at the delta function position we consider the infinitesimal interval x E [— e, e] obtaining 

0(e)=T.e/-^^(^)'^>(-e). (63) 

The integral is dominated by the delta function and for the problem we are treating in this 
paper we obtain the following boundary condition 

0(e) = e"'^"""^(/)(-e) . (64) 

To evaluate how the exponential acts on 0(— e) we use the Lagrange-Silvester formula PTll for 
a function f{M) of a matrix M 

;(M)^/(AO^ + /(A.)i^. (65) 



where Ai,2 are the eigenvalues of M. For the problem at hand, Eq. (|65l) leads to the following 
boundary condition around x = 



cosh(7 H-^sinhgf , (66) 






where g = is an adimensional interaction constant and 0"*^ represent positive and negative 
infinitesimals. A similar boundary condition holds for x = d. For the potential (l58l) we 
can now define three different regions, I, II, and III, defined asa;<0,0<x<(i, and 
X > d, respectively. In each of these regions the wave function is a sum of two plane 
waves of opposite momentum along the x— direction, with each plane wave multiplied by 
the coefficients Ar and Br, where T =1, II, and III labels the three regions defined above. 

Once the matrix T has been computed (see appendix [Appendix A[ ) the reflection 
coefficient is obtained from 

-^12 //—ix 
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Figure 12. (color on line) Top panels: transmitted flux Tf{E = 0.1, 0) for fixed g = 2 and 
two d =100a, 200a and for fixed d — 100a and three g =0.5, 1, 2. Bottom panels: transmitted 
flux Tf{E, 6 = 0) for the same previous cases. 



and the transmitted flux Tj from 

Tf = l-rr*. (68) 
For the case of zero electrostatic potentials, Ui = 0, we obtain 



^ 1 + sinh^ {2g) cos^ {2dk) ' 

which is a similar expression to the one in (|53l) if ^ = aw. It is simple to identify the limit 
T — ^ 1, which occurs when 2kd = {2n + l)7r, with n = 0, 1, 2, . . .. 

In Figure [121 we show the transmitted flux Tf{E, 6) as function of the energy and of the 
angle 6. The barrier between the leads and the center of the device is represented by a delta 
function potential, therefore wider barriers are represented by larger values of g. From Fig [12] 
we can see that for larger values of g the transmission in the forward direction is essentially 
zero except at some resonant energies, where the transmission goes to one. As function of 
the angle we see that there are some angles for which the transmission in also one. When the 
length of central part of the device is increase (larger d) the resonances become closer to each 
other and more resonances appear. 

In Figure [13] we present an intensity plot of Tf(E, 9) for a device with g = 0.5 and 
d = 200a. In this figure we can follow the evolution of the resonances in the E versus 9 
plane. The six lines of larger transmission are associated with the resonances we see in Fig. 
[l2]for g = 2mdd= 100a. 

As before, the linear response conductance per unit of cross-sectional length is given, at 
zero temperature, by Eq. (|55]) and represented in Fig[14] For small values of g the conductance 
shows smooth oscillations whereas for larger g values strong resonances are observed. The 
number of observed resonances depends on the length d. 
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Figure 14. Linear response conductance per unit of cross-sectional length at zero temperature. 
The top panels show G(/i, d) with g=0.5, 1 and the lower panels show the same quantity but 
forg=2, 4. 



8. Final remarks. 



In this paper we discussed the tunneling properties of Dirac electrons in two-dimensions 
when they transverse regions of space where the spectrum presents a finite energy gap. In 
the case we considered here the gap is induced by depositing graphene on top of Boron 
Nitride, rendering, in this way, the sub-lattices A and B non-equivalent. The consequence 
is the opening of a gap in the energy spectrum, that we have parametrized by the parameter 
t' = O.lt. We have shown that the existence of an energy gap prevents the Klein paradox from 
taking place, a necessary condition for building nanoelectronic devices made of graphene. We 
have also shown that basic devices like a resonant tunneling diode can be made of graphene, 
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by intercalating two regions where the spectrum of graphene presents a gap. We have also 
shown that simple analytical expressions can be derived for the tunneling through these types 
of heterostructures. In addition we have showed that a limiting case of the resonant tunneling 
diode can be understood by using Dirac delta function potentials. 

Clearly that one is lead to think that a full description of the ballistic (no impurities) 
transport process in the system should also include the effect of temperature and phonons. 
We note, however, that the electron-phonon interaction has been shown to have a small effect 
in the optical conductivity of graphene. [49J Which means that phonons should not be very 
important in the description of the transport process. Also the polaronic effect leads to a 
renormalization of the velocity vp vii the k — region and to the renormalization of the effective 
mass mt>|, = t' in the q — region. But since these two parameters can be considered as 
effective ones there is no point in including the polaronic effect explicitly. In concerns the 
temperature, clearly it will be of no importance when the chemical potential is above the gap. 
When the energy is in the gap there will certainly be temperature activated transport adding 
on top of the tunneling current. For small temperatures this will be a small effect. 

We believe our results of relevant for future nanoelectronics applications of graphene. 
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Appendix A. Details of section |7] results. 

We now give the details that allows one to derive the result [69] Applying boundary conditions 
derived in Sec. |7]and a T— matrix description ll42l for the scattering problem we obtain 




(A.l) 



with 



(A.2) 



The several matrices involved in Eq. (IA.2I) are defined as 





(A.3) 




(A.4) 
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+ tsinhg \^ ) , (A.5) 

and 



The momenta k2 and Ajs are given by 

1 



h = —f^\J{E - -{E- f/i)2sm^^i , (A.7) 

with i = 2, 3. The angles 9i are defined as 

k 

Oi = arctan -p , (A.8) 

with ky = \E — Ui\ sinOi/ [vfK)- The Sj functions are given by Si = sign(i? — Ui), with 
i = 1,2,3. The potential energies Ui represent some electrostatic potential created in the 
corresponding region. Although an analytical expression for T can be produced by carrying 
out the four matrix multiplications, the resulting expression is too cumbersome to be given 
here. In the special case that all Ui = 0, the matrix elements of T have a simple form given 
by 

Tn = cosh^ g + e-^''^'' sinh^ g , (A.9) 

and 

Ti2 = -sie-'^{l + e-'^'^'^) cosh ~g sinh g , (A.IO) 
with s = sign E, k = \E\ cos 6'/ {vph), and 9 the incident angle in the barrier. 
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